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In many architectures for fault tolerant quantum computing universality is achieved by a combination of 
Clifford group unitary operators and preparation of suitable nonstabilizer states, the so-called magic states. 
Universality is possible even for some fairly noisy nonstabilizer states, as distillation can convert many noisy 
copies into fewer purer magic states. Here we propose novel protocols that exploit multiple species of magic 
states in surprising ways. These protocols provide examples of previously unobserved phenomena that are 
analogous to catalysis and activation well known in entanglement theory. 
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Quantum computers are capable of executing algorithms 
whilst tolerating modest rates of faults or errors. Stabi- 
lizer codes encode information in subspaces of larger Hilbert 
spaces and allow a proportion of errors to be actively detected 
and corrected JT]- Whereas some anyonic systems with topo- 
logical^ protected ground states provide a passive method of 
safely storing quantum information (2)- Research into any- 
onic systems has been stimulated by the recent discovery of 
alloys that are topological insulators |3 , 4|, opening up a vari- 
ety of readily available systems that may be suitable for any- 
onic quantum computing. 

However, fault tolerant quantum computing is not just 
about archiving quantum information, but also processing the 
information whilst stored in its protected form. However, by 
employing stabilizer codes and topological systems we restrict 
how the quantum information may be fault-tolerantly manipu- 
lated. Stabilizer codes only allow coherent implementation of 
a limited group of fault tolerant gates, the so-called transver- 
sal gates. Unfortunately, recent research has shown that no 
stabilizer code can both protect against generic errors and of- 
fer a universal set of transversal gates [5|. Similarly, topo- 
logical^ protected groups of gates, implemented by braiding 
anyons, are not universal for many species of anyons JSHH). 
Theoretically, some exotic anyons do offer universal topolog- 
ically protected gates, but these are more physically specula- 
tive [9). Consequently, an alternative route to universal and 
fault tolerant quantum computing must be sought out. 

This obstacle is overcome by gate injection techniques. A 
suitable resource state is identified, and through fault toler- 
ant gates and measurements, this resource is consumed in ex- 
change for a new fault-tolerant unitary operator that promotes 
the group of gates to full universality. For both stabilizer codes 
and anyonic systems, the manifestly fault tolerant gates are of- 
ten contained within the Clifford group, the group of unitary 
operators that conjugate the Pauli operators. What resource 
states might promote the Clifford group to universality? Since 
the Clifford group maps stabilizer states — eigenstates of 
Pauli operators — to other stabilizer states, and such evolu- 
tions are efficiently classically simulable |10|, we know that 
stabilizer states fail to provide universality. However, numer- 
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ous nonstabilizer states do provide universality, including all 
single-qubit pure nonstabilizer states HI 1H . Bravyi and Kitaev 
proposed the appellation magic states for such resources fl2l . 
In their seminal article, Bravyi and Kitaev showed that some 
mixed nonstabilizer states can enable universal quantum com- 
puting via a process of distillation into purer magic states. 
Since preparation of the raw resources is not fault tolerant, 
we expect them to be noisy, and so distillation is essential. 

Some fault tolerance schemes actually provide a proper 
subgroup of the Clifford group, such as when braiding Ising 
anyons JSHH. Universality may still be possible via two lev- 
els of distillation if a resource state is available that first pro- 
motes the subgroup to the full Clifford group. For example, 
Bravyi [ 8 1 has shown that the aforementioned Ising anyon sys- 
tems can be promoted to the full Clifford group by distilling 
certain noisy stabilizer resources. 

The paradigm of magic states as a resource for promot- 
ing the Clifford group is analogous to other resource theo- 
ries, such as: how entanglement is a resource when only lo- 
cal operations are available 1131 : and how continuous vari- 
able Gaussian entangled states can be utilized provided with 
just local Gaussian operations fl4l . In both these alterna- 
tive examples of resource theories we have a thorough un- 
derstanding of the fundamental principles behind what state 
transformations are possible. The role of magic states is not 
yet understood as comprehensively as entanglement, although 
lately several results have begun to illuminate the subject. Re- 
ichardt ifTTl H31 [161 provided several additional distillation 
protocols beyond those found by Bravyi and Kitaev. He also 
identified some multi-qubit nonstabilizer states that can not, 
even probabilistically, be reduced to a single-qubit nonstabi- 
lizer state fTT|. Howard and van Dam IfTTl [T8l studied the 
role of noisy unitary operators as resources. They found that 
all depolarized single-qubit unitary operators that fall out- 
side the Clifford group can enable universal quantum com- 
puting. Campbell and Browne |[T9l l20l identified an analog 
to bound entanglement, with certain families of nonstabilizer 
states being undistillable for finite-sized computers. Ratanje 
and Virmani [21] considered resource theories that interpolate 
between separable states and stabilizer states and found new 
regimes that are efficiently classically simulable. 

This article explores the fundamental principles that govern 
magic states, and we uncover several new phenomena previ- 
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ously not observed. Many of the new phenomena have anal- 
ogous, though subtly distinct, counterparts in entanglement 
theory, such as entanglement catalysis 11221 and entanglement 
activation 11231 . Previous work on magic states has focused 
on what is achievable with many copies of the same quantum 
state. Whereas, a unifying theme of the protocols introduced 
here are the counterintuitive ways that two different sorts of 
resource can be jointly exploited. 

Magic catalysis can be described as a scenario involving 
two agents: a "magic state banker"; and an operator of a com- 
puter capable of only Clifford group operations. The banker 
is willing to loan magic states to the operator, but requires 
that the operator returns exactly the same quantum state at a 
later time. We identify a protocol where the loaned magic 
state acts as a catalyst, enabling the operator to perform state 
transformations that would have been impossible otherwise. 
Our protocol counteracts the misleading but intuitive idea that 
resources must be consumed to serve a function. 

Magic activation again involves a special resource, this time 
called the activator, that enables a probabilistic transformation 
that was impossible without this assistance. This phenomena 
differs from catalysis in several key ways. The activator is not 
returned to a banker, and the transformation may succeed with 
nonunit probability. Furthermore, the probabilistic transfor- 
mation also consumes a supply of bound magic states lfT9ll20l 
that alone have limited computational power when in finite 
quantity. 

Next we discuss the existence of, the aforementioned, com- 
putationally weak multi-qubit states that were first identified 
by Reichardt ifTTl . which we call irreducible non-stabilizer 
states. The defining feature of irreducible non-stabilizer states 
is that, on their own, no single-qubit nonstabilizer state can be 
extracted from one copy. We present new examples of irre- 
ducible non-stabilizer states for any number of qubits above 
two. Next we introduce another new protocol that exploits 
a combination of irreducible non-stabilizer states and bound 
magic states. Despite both resources being of limited util- 
ity we can, with some probability, extract a magic state of 
arbitrarily high fidelity. In many ways this protocol is more 
surprising than the previous magic state activation protocol. 
However, this latter protocol relies on a large number of re- 
sources. Depending on your preferred definition of activation, 
this protocol may also qualify as such. However we prefer 
to stress its unique aspects and so refer to it as an asymptotic 
activation protocol. 

Combined, these results provide a significant step towards 
a complete understanding of the principles governing magic 
states and their manipulation. Our results also prompt several 
interesting open problems that we discuss in the final section. 



I. TECHNICAL PREAMBLE 

In this section we refine our terminology and define nota- 
tion, beginning with a quick review of stabilizer states and the 
Clifford group. An rt-qubit pure stabilizer state, | ift) , is a quan- 
tum state uniquely defined by n commuting, and independent, 
Pauli operators gj . These operators generate by multiplication 



a group S of order 2™, the so-called stabilizer group for 
Every element of this group is said to stabilize the quantum 
state, such that s\ip) — \ip),\/s € S. More generally, a mixed 
state is a stabilizer state if and only if it is an incoherent mix- 
ture of pure stabilizer states. The Clifford group is the group 
of unitary operators that conjugate Pauli operators, such that 
for all Pauli operators p we have CpC^ = p' . Equivalently, 
the Clifford group are the unitary operators that preserve the 
set of pure stabilizer states. Important single-qubit Cifford 
unitary operators are the H (Hadamard) and T gates, which 
are best described in terms of their action on Pauli operators 

HXH^ = Z ; HZH^ = X; 
TXT^ = Y ; TYT^ = Z. (1) 

All single-qubit Clifford unitary operators can be decomposed 
into some sequence of these gates; that is, they generate the 
single-qubit Clifford group. To generate the entire multi-qubit 
Clifford group we have to add an entangling gate, such as the 
well known control-not gate. For further information on sta- 
bilizer states and the Clifford group, we refer the reader to 
Refs. mi24l. 

Throughout we refer to a Clifford computer as follows. 

Definition 1 A Clifford computer is a device capable of per- 
forming ideal Clifford unitary operators, preparation of sta- 
bilizer states, classical feedforward, classical randomness, 
Pauli measurements. 

For transformations implemented on such a device. 

Definition 2 If a Clifford computer can take an input state p 
and deterministically output a state p', then we denote this 
as p — >r> p', and say that p can be deterministically Clifford 
transformed to p' . Conversely, if there exists no such Clifford 
transform, we denote this as p p' . 

More generally, transformations may be probabilistic, as fol- 
lows. 

Definition 3 If a Clifford computer can take an input state p 
and with nonzero probability output a state p' , then we de- 
note this as p — >p p' , and say that p can be probabilistically 
Clifford transformed to p' . Conversely, if there exists no such 
probabilistic Clifford transform, we denote this as p -»p p' . 

The phenomena of catalysis and activation are essentially con- 
cerned with deterministic and probabilistic transformations 
respectively. 

The two most important single-qubit magic states are the 
eigenstates of the Clifford group unitary operators defined ear- 
lier, H and T, such that 

H\H ) = \H ) ; H\Hx) = -\Hi); (2) 
T\T Q ) = e i7r / 3 \T ) ; T\Ti) = e _i,r / 3 |Ti). 

We also use similar notation for stabilizer states such as Y 
eigenstates | Yo,i) ■ For an n-qubit state a binary vector v = 
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{vi, ■■■v n ] specifies the state 



\H v )=(g)\H Vj ), 



(3) 



and similarly for |T V ). Employing greek characters for mixed 
density matrices, we use 

t v = \T v )(T v \, (4) 

with v again an n-bit vector. 

II. MAGIC CATALYSIS 

Here we present an example of magic catalysis. 

Theorem 1 Magic catalysis is possible: for the state \cp) oc 
|-ffo,o,o) + we have \cp) \Hq) but with the addi- 

tion of catalyst \Hq) we have \<p)\Ho) — >d \Hq)\Hq). 

Clearly, this satisfies the constraints of the scenario described 
in the introduction since the process is deterministic and the 
catalyst is unchanged it can always be returned to the banker. 
First we describe a protocol, also illustrated in Fig. [T] that 
implements the deterministic transformation \ip)\Ho) —>d 
\H )\H Q ). 

1. Prepare the state \<p) on qubits A, B, C, and state \Hq) 
on qubit D; 

2. Measure the Pauli stabilizer YcYd; 

3. If the measurement yields outcome +1, then apply the 
unitary operator Hp; 

4. Measure the Pauli stabilizer ZqZd; 

5. If the previous measure yields outcome — 1, then apply 
the unitary operator YaYb; 

6. keep qubits A and B, and discard qubits C and D. 

Although the process involves two measurements with ran- 
dom outcomes, each measurement is conditionally followed 
by a unitary operator that ensures the same output regardless 
of the measurement outcome. Consider step 3, after the Yc^d 
measurement with a +1 outcome, we have the state 

H D (i+Y c Y D )\<p)\H ) = (l-Y c Y D )\<p)H D \H ), 

= {t-Y c Y D )\cp)\H Q ), (5) 

where the first line uses HjYj — —YjHj, and the second line 
uses H\Hq) = \Hq). Hence, we can deterministically imple- 
ment a projection onto the —YcYn subspace. 

Next, measurement results —ZqZjj or -\-ZcZjy give a pro- 
jection of these qubits onto the state \^^) oc |1, 0) — |0, 1) or 
|$ + ) oc |0,0) + |1, 1), respectively. The use of projec- 
tions plays a pivotal role throughout this article, effectively 
functioning as an odd parity projector for any basis. That 
is, for any orthonormal basis {|&o)> IMI shared between two 
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FIG. 1: An outline of the magic catalysis protocol. Filled circles rep- 
resent qubits, and lines between qubit denote correlations. The pro- 
tocol involves two measurements with random outcomes. Although 
different measurement outcomes produce different projections, adap- 
tively applied Clifford unitaries ensure the outcome is deterministic. 
However, the determinism of our protocol relies heavily on the sym- 
metries of the initial states. The quantum states \tp) and \Ho) are 
defined in Thm.[T]and Eqn. {2|, respectively. 



qubits we have = (1 — 5j^)/v2, where 5j t k 

is the Rronecker delta. This feature of the singlet projector 
follows from (U ® U)\^>^) oc for any unitary opera- 

tor U, and so l^ - ) is odd-parity in any basis. Returning to 
the problem at hand, the relevant basis is the Hadamard basis, 
where \^f~) oc |i?o,i) — l-ffi.o)- Hence, the singlet projection 
picks out the second term of |i?o,o.o,o) + |-ffi,i,i,o)> producing 
\^~), In accordance with step 5, we apply YaXb (not- 
ing F-gates flip Hadamard eigenstates) and discard the last 
two qubits. This yields the desired output |i?o o)- 

If instead, step 4 provides a -\-ZcZd measurement out- 
come, we have a projection onto the state |$ + ), and so 

($+\ CtD \<p)\H ) ex {*-\c,DY D {\H m o) + |ffi,i,i, », 
oc (*-|c,.d(|#o,o,o,i) + 1-^1,1,1,1)). 
<* |i?o,o), 

where the first line uses |$ + ) oc Yd\^~), allowing further 
employment of the singlet projector. Hence we yield the de- 
sired output regardless of measurement outcomes. 

To prove that we have identified a truly catalytic process, 
we must also show that the process was otherwise impossible, 
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such that \<p) -i*d \Hq). We actually proceed by showing the 
stronger result that \ip) -»p \Hq), which directly entails the 
weaker deterministic no-go result. Since we are attempting to 
probabilistically output a single-qubit state, we only have to 
consider Clifford transformations that project onto a stabilizer 
codespace, with a single logical qubit, and then decode lfl9l . 
For a codespace with logical states |0l) and the result of 
projecting and decoding performs the transformation 



|y> ^ IVW) cx <0i|^>|0> + <li|y>)|l) 



(6) 



For projections onto stabilizer subspaces, the ratio of the com- 
putational amplitudes, 



fl(IVw t » = KQiMlVKUMI 2 , 



(7) 



must be one of a few possible rational fractions (see Ap- 
pendix [A}. However, for the target state, \Hq), this ratio is 
an irrational number tan 2 (7r/8) = 3 — 2y/2. Hence, the exact 
transformation is impossible, and our proof is complete. 

The techniques in Appendix [A] are sufficiently general to 
rule out many other Clifford transformations. For example, if 
we ask whether n copies of \tp) can be exactly converted into a 
instate, our method also proves this is impossible with finite 
n, and so \ijj)® n -Prp \Hq) (discussed further in Appendix|B|>. 

In the preceding example of catalysis, both initial quantum 
states are pure. However, there is a variant of the protocol 
where \ip) is replaced by a mixed state, 



(\<p)(<p\+H**\<p)(<p\H 



(8) 



It is easy to prove that catalysis can be performed with a v be- 
cause it is Clifford-equivalent to \(p); that is |<^)(<y5| -^d °V 
and <r v — >d I^X^I- The first transformation is clearly 
possible by randomly applying H® 3 to \ip). The converse 
transform is achieved by measuring the observable YaYbYc- 
Given the +1 outcome we have the state \(p), whereas the —1 
outcome produces H® 3 \ip). 

Describing catalysis in terms of an interaction between a 
magic state banker and a computer operator gives it an op- 
erational flavor. Although the scenario could be considered 
somewhat artificial. We feel that the true depth of catalysis 
is that certain transformations become possible, for free, as- 
suming a reserve of magic-states. Such transformations can 
be called magic-assisted Clifford group operations, and it is 
an interesting open problem to determine the full structure of 
these operations. 



in. MAGIC ACTIVATION 

Here we give an example of magic activation. One of the 
distinguishing features of activation is that it utilizes resources 
from a family of bound states. Rather than a general account 
magic-state boundness, for brevity we describe the concept 
with respect to noisy T-magic states, 
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FIG. 2: An outline of the magic activation protocol which exploits 
the bound state r(/) (see Eqn. 9} in the presence of its activator a(q) 
(defined in Thm.[2|. This single shot protocol succeeds probabilisti- 
cally when qubits B and C are projected onto the singlet state. 



where / is the fidelity and / st = (1 + l/\/3)/2 is the thresh- 
old above which we have a nonstabilizer state. The following 
statement follows directly from the more general results of 
Ref. 



Theorem 2 For any finite n, there exists a positive e n > 0, 
and a corresponding no-go region of fidelities f < f st + en- 
Inside this no-go region, it follows that for any single qubit 
state, p, we have that r(/)® n — >p p if and only if t(/) — >p p. 
We say that the family of states r(/) is bound. 

Heuristically, this result instructs us that there exist nonstabi- 
lizer states where n copies are no more useful than a single 
copy. Since this holds even with probabilistic postselection, 
we cannot distill these states to higher purity. There is clearly 
a parallel with bound entanglement, but there is also a subtle 
distinction. The threshold fidelity, / st + e„, below which the 
theorem applies, depends on the number of copies, n. Hence, 
it is possible that the region shrinks as n is increased, maybe 
even such that e„ — > as n — > oo. In contrast, bound en- 
tangled states are bound regardless of how many copies we 
have. However, it is not known that the region actually does 
shrink. Rather, it is merely a limitation of the techniques of 
Ref. 11201 that this possibility has not been ruled out. Whilst 
known techniques |[T2l can distill noisy T-states with fidelities 
greater than (1 + W3/7)/2, there is no known method of dis- 
tillation that functions below this fidelity. Hence, it is possible 
that even for large n the no-go region does not shrink below 
this level. 

Subtleties aside, it is clear that r(/)®", with sufficiently 
small fidelity and fixed n, cannot be distilled. This is in con- 
trast with noisy instates, which are not a bound family of 
states. For example, consider noisy instates with any initial 
fidelity large enough that no stabilizer decomposition exists. 
With 7 copies of such noisy instates one can implement a 
protocol lfTTI[T6l based on the STEANE code that, when suc- 
cessful, increases the fidelity ll27ll . The protocol must be iter- 
ated to achieve higher fidelities, and a unit fidelity is asymptot- 
ically approached with increasing n. However, the important 
feature is that some fidelity increase is always possible with 
finite copies, and that a similar protocol for all noisy T-states 
is ruled out by Thm. [2] This prompts the question, are very 
noisy T-states ever useful resources! We affirmatively answer 
this question by providing an activation protocol. 



r(/) = /r + (l-/)n, 



(9) 
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Theorem 3 Magic activation is possible: for the activator 
a {<i) — Q T o,i + (1 — <l) T i.o (f or some 1 > q > 1/2) and 
any r(/), with f st < /, there exist a single-qubit state p such 
that: 

i. o~(q) ® t(/) — >p p; even though 

ii. <r(q) p; and 

Hi. t(J) ^> p p. 

Alone neither state can produce a particular output p, but com- 
bined it is possible. The output state is again a noisy T-magic 
state, so p = r(f'). Provided /' > /, condition (2'2'z) of the 
theorem immediately follows. 

We begin by describing the activation protocol, also illus- 
trated in Fig. [2] 

1. prepare state a(q) on qubits A, B and r(/) on qubit C; 

2. measure the observables YbYq and ZpZc', 

3. postselect on —1 for both measurement outcomes, and 
discard qubits B and C. 

The initial state can be expanded out as 

<q)®r{f) = (7/70,1,0 + (1 - g)/ri,o,o (10) 

+ ?(i-/Ki,i + (i-g)(i-/h,o,i. 

The postselected measurements project qubits B and C onto 
the singlet state. We use that |^~) oc |Ti, ) - |T .i), and so 

(*>(</) ®r(/)l*~) cx qfr + (1 - <?)(1 - /)n(ll) 

We have effectively projected onto the odd parity terms of 
qubits B and C, and then traced them out. After renormal- 
ization, the state is r(/') with fidelity 

f> = tl (12) 

1 q..f + (l-q)(l-f) ^ ; 

It is easy to see that /' > / whenever 1 > q > 1/2, and 
so the transformation could not be achieved with r(/) alone, 
satisfying condition (Hi). To complete the proof we must show 
condition (ii); that the transform could not be achieved with 
a(q) alone. 

The simplest transformation on a(q) alone is to measure 
qubit A of a(q) in the computational basis. Due to the T- 
symmetry of the state, any Pauli basis gives the same result. 
Hence for a ±1 outcome of any single qubit Pauli measure- 
ment, the resulting unnormalized state is 

ti A [(1 ± Z A )a(q)] cx qc±n + (1 - - c±)t , (13) 

where 

c± = trace [(1 ± Z a )t ] /2 = (l ± /2. (14) 

Clearly, the outcome gives a greater fidelity. Further- 

more, we have c + = / st . Renormalizing gives a noisy T-state 
with fidelity 

f" = ^ (15) 



This fidelity fails to match that achieved by our activation 
protocol. However, a single qubit observable is clearly not 
the only option available, with many possible stabilizer mea- 
surements over both qubits. Checking other possible mea- 
surements (see Appendix [C| one finds that the simple single- 
qubit measurement proves to be optimal. Therefore, a single 
copy of a(q) cannot be probabilistically Clifford transformed 
to p(f'), the output of the protocol, and so the activation is 
genuine. Of course, our argument does not rule out that many 
copies of a(q) may accomplish this transformation, as is in- 
deed the case (see Appendix|E|. This feature is consistent with 
the analogous phenomena of activation in entanglement the- 
ory 1123 1 , as known entanglement activators are also many- 
copy distillable. 

It is unclear whether more copies of the bound resource 
could be exploited to iterate or improve this particular magic 
state activation protocol. However, the subsequent sections 
describe a more involved protocol that is stronger in two prin- 
ciple respects: firstly it can be extended to consume arbitrarily 
many bound resources, with an output fidelity asymptotically 
approaching unity; secondly, the activating resources is also a 
computational weak state of a species that we introduce next. 

IV. IRREDUCIBLE NON-STABILIZER STATES 

This section introduces the notion of an irreducible non- 
stabilizer state, which is another form of noisy resource that is 
computationally weak. We also present some new examples 
of such states to be used in the next section. 

Definition 4 A state a is an irreducible non-stabilizer state 

(an INS state) if both 

1. a is not a stabilizer state; and 

2. for all single-qubit nonstabilizer states, p, we have 
a p. 

Reichardt identified the first examples of INS states ATI . 
However, Reichardt referred to them as counterexample 
states, as he presented them to disprove a conjecture that all 
multi-qubit nonstabilizer states can be Clifford transformed 
to a single-qubit nonstabilizer state. Obviously there are no 
single-qubit INS states, but Reichardt showed that two-qubit 
INS states do exist. Despite being of limited computational 
power, some INS states prove useful when combined with 
other resources. We consider states of the form 

a ws (q,n) = qT® n + (l-q)t/2 n ; (16) 

which satisfy the definition of an INS state whenever the 
weighting, q, falls in a specific interval, q min < q < q max , 
where 

<Zmax = [l + (2/ st )"- 1 (y3-l)]- 1 , (17) 

g min = (2"-l)/[(l + x/3)"-l]. (18) 
Values of q and n satisfying these conditions are shown in 
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FIG. 3: A region of INS states of the form ctins (q, n ), as in Eqn. l )16} . 
The weighting, is shown on a log-scale, against the number of 
qubits n. All states with the weighting satisfying q m in < q < g max 
are INS states, with the region being empty for n — 1,2 and ap- 
pearing for n > 3. Some states outside the shaded region may also 
qualify as INS states. 



First we show that for sufficiently pure states, we indeed 
have a nonstabilizer states, and so meet condition (1) of the 
definition. In general, mapping out the space of multi-qubit 
mixed stabilizer states is an involved problem ifTTll . However, 
there exists a simple witness that can detect many nonstabi- 
lizer states. We introduce this witness in terms of a norm we 
call the stabilizer-norm (or just st-norm): 

Lemma 1 A density matrix p, with decomposition in the Pauli 
basis p — djUj, is a nonstabilizer state if 
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a,- > 1. 



(19) 



For single qubit states the condition is not just sufficient, but 
also necessary. Indeed, for a single qubit state this inequality 
marks out an octahderon in the Bloch sphere. However, there 
are many multi-qubit nonstabilizer states that are not detected 
by this witness. To prove the lemma we first observe that the 
st-norm satisfies the triangle inequality, and hence is convex. 
Furthermore, all pure stabilizer states, p st , have unit st-norm, 
||/?st||st = 1, and so no mixed stabilizer states can exceed 
unity. 

For the states of interest here, the st-norm is 



lkiNs(9,n)||s 
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+ (1 - <7)/2 r 
(l-?)/2", 



(20) 
(21) 



where the second line uses multiplicity of the st-norm with 
respect to the tensor product; in general \\p a g3 pbllst = 
I I Ut- 1 1 1st- Calculating ||r || s t = (1 + VS)/2, and re- 
quiring the st-norm exceed unity, entails q > q m i n . 

Next we prove that for sufficiently impure states, condition 
(2) of our definition is satisfied. It is well known |[T9l that 
such a transformation is impossible if it cannot be achieved by 
projecting onto a stabilizer codespace, with one logical qubit, 
and decoding. First we note that all mixed single-qubit states 



with largest eigenvalue satisfying A < / s t are stabilizer states. 
We prove that, for q < q niax , all codespace projections fail to 
achieve sufficient purity. Hence, they output stabilizer states. 
For a stabilizer code with projector n, the projected state is 



Pout — 



g.nrf "n + (i - g )n/2" 

g.tr(nr ®") + (1 - 9 )/2«- 
The largest eigenvalue of the projected state is 

g.tr(nr ®") + (1 - g)/2» 



A 



q.tv(IlT® n ) + (1 - q)/2»- 



(22) 



(23) 



To make further progress we must evaluate the maximum pos- 
sible value of tr(nr ®"). 

Lemma 2 For n copies of a single-qubit state, tq, and for all 
projectors, II, onto a 2 m -dimensional stabilizer subspace, the 
maximum probability of projection is 



max[tr(nrr)] = /, 



n—m 
st 



(24) 



This lemma asserts that the maximum probability of any sta- 
bilizer projection is achieved by a series of single qubit sta- 
bilizer measurements. The lemma can be proven using graph 
codes [25] as shown in appendix|D] Applying the lemma (with 
m = 1) to Eq. (23 i gives a maximum achievable value of A, 
which we denote with a star 



A* 



qf^ + jl-q)/^ 
qfST 1 + (1 - q)/2*-i 



(25) 



If we wish to guarantee that the output is a stabilizer state, 
we require A* < / s t, and a little rearrangement produces the 
inequality q < q max . 

Hence, we have proven the existence of a whole class of 
INS states using a very different approach to Reichardt. Note 
that q m i n and q max differ only for three or more qubits, so our 
construction does not provide any two-qubit INS states. 



V. ASYMPTOTIC MAGIC ACTIVATION 

Another new protocol is described here, which demon- 
strates two features not exhibited by the previous activation 
protocol. Firstly, its exploits a combination of an INS state 
and many bound magic states. Secondly, the output magic 
state can be arbitrarily pure, asymptotically approaching unit 
fidelity with the number of bound states used. In light of this, 
we distinguish this protocol by calling it asymptotic activa- 
tion. For the purposes of this section, we consider the INS 
states with q = g max , as in Eq. ( 20 1, and for brevity herein 
use the notation 



OINSM = CMS (^max,™)- 

Using this resource we have the following result. 



(26) 



Theorem 4 Asymptotic magic activation is possible: for the 
INS state ctinsM and any t(/) (8 ™~ 1 (with f > f st ), we have 
that cr IN s( n ) ® t(/)' 8i ™ _:L ->p T~(f), where /'->lfljn-> 
oo. 
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FIG. 4: An outline of the asymptotic magic activation protocol. The 
protocol uses n — 1 copies of the noisy T-states r(f) (see Eqn.[9]( and 
a specific n-qubit activator ams(q) (see Eqn. [26]). The activator is 
an especially weak resource, known as an irreducible non-stabilizer 
state (defined in Def. |4). The protocol pairs up each noisy T-state 
with a qubit from the activator, and succeeds when all n — 1 pairs are 
projected onto the singlet state. 



By definition, the INS state cannot be reduced to a single- 
qubit nonstabilizer state. Since the bound states also resist dis- 
tillation, the protocol seems to exploit some synergy between 
the two resources. As with the previous activation protocol we 
utilize singlet projection. The asymptotic activation protocol, 
also illustrated in Fig.|4j is as follows. 



1. prepare <7iNs( n ) on qubits A,B, 
T{f)® n - 1 onqubits A',B',...; 



and prepare 



2. flip every qubit of r(/)® n 1 using the local Clifford 
H.Y that maps tq — > t\\ 

3. pair up n — 1 qubits from each resource, pairing A with 
A' and B with B' , etc; 

4. Measure the observables XjXj> and ZjZf for every 
pair; 

5. Postselect on outcomes for every measurement 
outcome in every pair, and discard all measured qubits. 



After step 2, the quantum state is 



5n-l 



P = gmaxTfP 65 r(l - ff 

-(l-gm»x)^®r(l-/) 5 



(27) 



in — 1 



The subsequent steps project on the singlet state between 
paired up qubits, giving 



P oc g max a 



L ro + (1 



where 



a = (¥-|7rj®r(l -/)|*-> =//2, 
b = (*-|l®r(l-/)|* _ )/2 = l/4. 



(28) 



(29) 



Combining these equations and after some manipulation, we 
find that the output fidelity is 



/' = 



fat 

T 



(V3-1) 



(30) 



Clearly this approaches unity in the large n limit, provided 
that / exceeds the stabilizer threshold. 

Unlike the previous activation protocol we are allowing for 
the number of copies to vary, with the phenomena becom- 
ing more pronounced in the large n limit. Since the number 
of copies is varying, and not fixed to some finite n, previous 
no-go results on the non-distillability of t(/)®" do not ap- 
ply. Consequently, we cannot guarantee that the transforma- 
tion would be impossible without the addition of the INS state. 
As we have not strictly proven r(/) (8n -^p r(f'), we have 
exercised caution and not described this as a vanilla activation 
protocol. However, for small fidelities, / < (1 + y / 3/7)/2, 
there is no known protocol [12] that performs this transforma- 
tion, even in the limit of many copies. The lesson this protocol 
teaches us is that large numbers of noisy T-states can be ex- 
ploited to great effect when accompanied by another resource 
state. The most fundamental open problem in this research 
area is now whether asymptotically many nonstabilizer states 
can be purified when completely unassisted by activating re- 
sources. 



VI. DISCUSSION AND CONCLUSIONS 

We have introduced three protocols for quantum computers 
with Clifford-group unitary operators that are fault tolerant, 
and for clarity are taken to be ideal. All our protocols make 
use of two different species of nonstabilizer states, which is 
a relatively unstudied topic compared with that of distilling 
many identical copies of a quantum state. Each of the pro- 
tocols is designed to illustrate a peculiar and counterintuitive 
phenomena that can occur in Clifford computers. We now 
review each of these protocols and discuss related open prob- 
lems. 

Magic catalysis demonstrates that reserves of magic states 
do not have to be depleted to serve a function. A magic state 
can act as a catalyst that enables a deterministic transforma- 
tion that is impossible by Clifford transformations alone. In 
our catalytic protocol, the catalyst was a Hadamard eigen- 
state, and the protocol depended on some very specific sym- 
metries of this state. This prompts the question of whether 
other nonstabilizer states can serve as catalysts. For example, 
can eigenstates of the T-gate also act as catalysts? We conjec- 
ture that — in light of deep underlying differences between H 
and T-states — the answer will be no. The T magic states are 
weaker in several regards. Firstly, existing proposals for im- 
plementing non-Clifford gates do not directly exploit T-states. 
Rather the T-states are probabilistically converted into states 
on the Bloch sphere equator (see Ref. | fT2) or Appendix [Q , 
and only then are they used for implementing a non-Clifford 
gate. Secondly, noisy T-states just outside the set of stabilizer 
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states are undistillable, or bound, in the sense reviewed ear- 
lier. Beyond this anecdotal evidence, we have no firm proof 
that T-states cannot function as catalysts. However, settling 
the conjecture either way should prove illuminating. 

No protocol, prior to this article, has exploited noisy In- 
states arbitrarily close to the set of stabilizer states. Indeed, 
the evidence surveyed in the previous paragraph suggests that 
there exist noisy T-states, which despite being nonstabilizer 
states, cannot be utilized for any useful task. However, our 
magic activation protocol shows that a noisy T-state com- 
bined with an activator resource can probabilistically output 
a single-qubit state that could not be achieved with either re- 
source alone. Hence, all noisy T-states outside the set of stabi- 
lizer state are useful for some task. Since all noisy H nonsta- 
bilizer states are already known to be useful without the assis- 
tance of an activator, it is somewhat redundant to ask whether 
activation could be performed with H states |28|. A dissim- 
ilarity with entanglement activation is that our magic activa- 
tion protocol is not iterative, being defined for only a single 
round. The most interesting questions on this topic concern 
what kinds of iteration are possible. Our third and final proto- 
col, asymptotic activation, gives one possible extension. 

Asymptotic activation shows that (n — 1) copies of any 
noisy T-state and a particular n-qubit resource can probabilis- 
tically output a magic state, which in the asymptotic limit ap- 
proaches unit fidelity. Furthermore, the special n-qubit re- 
source is an irreducible non-stabilizer state, from which no 
single-qubit nonstabilizer states can be probabilistically ex- 
tracted. The class of irreducible non-stabilizer states is in- 
teresting in it own right, so our methods for constructing them 
may find applications elsewhere. Indeed, one interesting prob- 
lem is whether many copies of irreducible non-stabilizer states 
are distillable or a new form of bound state. 

Neither asymptotic activation nor standard activation are 
analogous to entanglement activation in every respect. For ex- 
ample, in the entanglement activation of Ref. [23 ] the protocol 
simultaneously exhibits the following three features: 

1. The protocol can consume a variable number, n — 1, 
copies of the bound resource with the output fidelity 
tending towards unity with increasing n; 

2. The activating resource has a fixed size; 

3. It is proven that neither the bound resources nor the ac- 
tivating resource can on their own be probabilistically 
transformed to the output of the activation protocol. 

Asymptotic activation has property (1), standard activation 
satisfies (2) and (3), but neither magic protocols simultane- 
ously exhibit all three features. This prompts the question, do 
magic protocols exist that are more sturdy analogs of entan- 
glement activation with all three features! The extent of sym- 
metries between the two resource theories is far from clear, 
and hence so is the answer to our posited problem. 

Considering all our protocols together, a key tool in all is 
the use of a singlet projection along with at least one state with 
multi-qubit correlations. The singlet projector functions as a 
method of verifying if two qubits are nonidentical, although 
at the price of projecting those qubits into a stabilizer state. 



Since our aim is to prepare nonstabilizer states, singlet projec- 
tions can only be exploited when accompanied by multi-qubit 
correlations. Indeed, we have seen that singlet projection is 
an extremely useful tool in this context. So far we have not 
considered any scenarios with many copies of a multi-qubit 
correlated state, but it seems plausible that the singlet projec- 
tion would prove useful in such contexts. This is indeed the 
case, and for completeness we provide just such a strategy in 
Appendix [E] 
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Appendix A 

Here we show that the ratio of amplitudes in Eq (|7]) can only 
take rational values, and hence cannot achieve the required ir- 
rational number. Furthermore, the set of possible ratios is fi- 
nite, so there is a limit to how closely the target ratio can be 
approximated. We present a very general form of the proof, 
which can be used to rule out many other Clifford transforma- 
tions. The techniques introduced here indicate that there may 
well be hope for building a general framework for understand- 
ing catalysis. 

Before beginning the core proof, we shall make some ob- 
servations; The specific initial state of interest, \ip), is Clifford 
equivalent to 

y) = (|0, 0, 0) + i|0, 1, 1> + 0, 1) + ijl, 1, 0))/2. (Al) 

The required Clifford is simply \/X , which maps the 
Hadamard eigenstates to the Bloch sphere equator such that: 

VX\Hi) = \H$) = (|0) + {-iye-^\l))/V2. (A2) 

Hence, an equal superposition of \H' Q o) an d \ ^\ l l) cancels 
out odd excitation terms, leaving only the even terms shown 
in Eqn. (| A1| >. Next we recall that, stabilizer states must have 
the form 11261 

KU= E * i(x) (-i) 9(x) |x + y)/v^ (A3) 

xGC+y 

where Z(x) and q(x) are some linear and quadratic functions 
in x, C is a binary linear subspace and y is some constant 
binary vector. Notice that our state \<p') has a very similar 
form to stabilizer states as in the computational basis the co- 
efficient have equal magnitude and phases are multiplies of i. 
Such states are interesting and deserving of their own title. 
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Definition 5 We say a pure quantum state, is a pseudo- 
stabilizer state if and only if there exists a Clifford unitary U 
such that: 



(A4) 



t/|<A) = (E j/(x) l x »/v^ 



where V is a set of n-qubit bit strings, and f : x — > {0, 1}. 
Furthermore, we say a pseudo-stabilizer state has complexity 
P, such that 



P{\ip)) = min ■ 



f _ 7 -/( x )| x 



■yu g c 



(A5) 

This is simply the smallest possible \V\ over all valid decom- 
positions. 

Notice that genuine stabilizer states also satisfy this definition 
but have trivial complexity P = 1, and the decomposition of 
Eqn. (|A1| entails that P(\<p')) < 4. In contrast the instates 



are not pseudo-stabilizer states as defined above. 

Here we prove that the complexity of a pseudo-stabilizer 
state limits the possible single-qubit states one can produce 
by Clifford transformations. 

Theorem 5 Consider a pseudo-stabilizer state of com- 
plexity P. If \tp) —}p \ipaut) where \ipout) is a pure single 
qubit state, then it follows that the amplitude ratio satisfies 



i?(IVw)) 



where Rp is the set of feasible ratios 
,2 



Ri 



'0 



aj,bj G Z; \ aj \ + \bj\ < P 



(A6) 



(A7) 



Conversely, for any \t])') with (fi Rp we can conclude 

\ip)^ P W). 

Notice that all of the feasible ratios from a pseudo-stabilizer 
state are rational fractions, so exactly producing a instate is 
impossible. The specific result \ip) -»p \Hq) then follows 
from our earlier observation that \ip) is a pseudo-stabilizer 
state of bounded complexity P(\<p)) < 4. 

As noted in the main text, we only have to consider prob- 
abilistic Clifford transforms that project onto a single qubit 
stabilizer subspace and then decode, and so we can achieve 



R(H out )) = \(o L m 2 /\(i L m 2 , 



(A8) 



where |0l) and are logical states of the stabilizer sub- 



space, which by Eqn. ( A3 1 can be expressed as 

j'o.l(x)(_ -|^go,l ( x )| x } 



0,l£ = 



xeC ,i+yo,i 



(A9) 



with the numeric subscripts differentiating (C, y, q, I) for the 
two states. It is well known that logical states of stabilizer 



codes can always be found such that they differ by Pauli rota- 
tions, such that \1l) = Xl\0l)- Pauli operators can change 
y, I and q, but not C and so Cq = C\ = C. 



Using equation ( A4 1 and ( A9 1 we find that: 



<o,w)= E 

xe-pn(c+y ,i) 



j/ ,i(x)+/(x) ^_2)9o,i( x ) 



(A10) 



each term in the summation is a multiple of i and there are no 
more than P(\tp')) — \V\ terms. Hence we have 



(0, l L |V/> = Ki + ibo,i)/y/\C\.\V\, 



(All) 



where aj ,bj G Z and the limited number of terms entails 
\a>j\ + \bj\ < P(\i/j')). Calculating the ratio of these am- 
plitudes, the ICI-lPl factors cancel and we have the result as 
stated in Thm. 

From an infinite set of rational numbers, one can always 
find an arbitrarily good approximation to any irrational. How- 
ever, the set of feasible ratios is limited by the constraints 
\ a j\ + \bj\ < P(\ip)) and and so the set of possibilities is 
not just finite but potentially very small. We have presented 
an argument based on rationality for generality. However, it 
is quite straightforward to numerically search the limited set 
of possibilities and verify that for P(\f)) = 4 we can never 
achieve R = tan(7r/8) 2 . Such a search produces 1/5 as the 
closest possibility, which differs from the target by over 0.028. 
Note that Thm. |5]l places a restriction on feasible ratios, but 
does not guarantee that all such ratios are achievable. 

The theorem deduced here rules out many Clifford trans- 
forms, but is far from the generality of the majorization crite- 
ria that is used in entanglement theory ll22l . In entanglement 
theory, the majorization criteria depend on the coefficients of 
the quantum state in the Schmidt basis. If we consider all 
possible local unitaries, and rotate a state to have the mini- 
mal possible support in the computational basis then this also 
yields the all important Schmidt coefficients. Returning to the 
context of magic states, our approach hints that minimizing 
support over all possible Clifford unitaries also gives a de- 
composition with important coefficients. Our investigations 
into this approach are ongoing. 



Appendix B 

Here we briefly address the question of whether many 
copies of \(p) can be probabilistically Clifford transformed 
into \Hq). Much of the technical apparatus required was es- 
tablished in Appendix [A| Given that \tp) is a pseudo-stabilizer 
state with complexity P[\ if)) < 4, it follows that | tp) ® n is also 
a pseudo-stabilizer state but with P(\ijj)) < 4". Hence, for fi- 
nite n, Thm|5]applies and we can conclude -Prp |i?o)- 
However, a supply of \(p) states is a resource for universal 
quantum computation. This paradox is resolved by observ- 
ing that although an exact \Hq) is impossible to produce, we 
may approximate \Hq) with a fidelity that asymptotically ap- 
proaches unity as n increases. Similarly, in entanglement the- 
ory many copies of any pure entangled state may be converted 
into any other state in the asymptotic limit. 
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Appendix C 

We consider two qubit stabilizer measurements on the state 
a(q), defined in Thm. [3] which is an incoherent mixture of 
|7o,i) an d |Ti,o) • We shall exploit that the state is invariant 
under T rotations of either qubit A or B, such that 

T%T b B a{q){T%T b B y = a(q) (CI) 

for any integers a and b. Furthermore, for any Pauli opera- 
tor P a Pb, where P AjB = {X A>B ,Y A<B , Z AB \, there exists 
integer a and b such that: 

T%T b B P A P B {T%T b B y = Z A Z B . (C2) 

Combing these two properties of the T rotation, we have that 
for any two-qubit Pauli projection 

T%7%(1 ± P A P B )a(q)(i ± P A P B )(T%T b B y 

= (i±Z A Z B )a(q){i±Z A Z B ). ^> 

This symmetry entails that we only have to consider two pos- 
sible Pauli projections, II± = (1 ± Z A Z B )/2. As an inter- 
mediate step in our proof, we see that when the state is pure, 
q = 1, two T-states can be probabilistically converted into a 
pure state on the Bloch sphere equator. This equatorization 
is an important step in using these resource for implementing 
non-Clifford gates. 

First, we note that T-states in the computational basis are 

|T ) = cos(/3)|0)+e M / 4 S in(/3)|l), 

|Ti) = sin(/3)|0)- e ^ 4 cos(/3)|l), (C4) 

where cos(2/3) = 1/V3. If we consider the n + projection 
onto the even parity subspace then 

n+|r 0) i) =n + |T 1)0 ) = cos( ( S)sm( ) (3)(|0,0)-i|l,l))- 

(C5) 

Since either pure state is projected onto the same stabilizer 
state, so too is the mixture <r(q). 

For the odd parity projector, n_, the analysis is more in- 
volved as 

n_|T 0ll > = e " /4 [sin 2 (/3)|l,0)-cos 2 (/3)|0,l)] , 
n_|T x ,o) = e " /4 [sin 2 (/3)|0,l)-cos 2 (/3)|l,0)] , 

which are distinct nonstabilizer states. Using the decoding 
|0, 1) — > |— ) and |1, 0) — > — these states map to points 
on Bloch sphere equator, 

n_|T 0ll > | 7+ ) = (|0) + e^|l))/A 
n_|T li0 > | 7 -) = (|0) + e-*T|l))/v^, 

where 7 = ir/6. Since I7+) is in the positive octant of the 
Bloch sphere, no other decoding gets closer to the target \Tq) 
state. Applying this analysis to the projection of the initial 
mixed state gives 

n_a(q)n_ <z| 7+ }(7+l + (1 - ?)|7->(7-|. (C6) 



The fidelity of this output with respect to |T ) is 

/"' = ?l(rol7 + )l 2 + (i-?)l(r |7->l 2 , (C7) 

where, 

l(To|7±)| 2 = (9±V3)/12. (C8) 

Comparing the fidelity /"' with /" of Eq. $15) , we find that 
/"' is always smaller. Hence, no two qubit stabilizer projec- 
tions can outperform the single qubit projection. 

Appendix D 

This appendix provides a proof of lemma [2] which gives 
the maximum probability of projection onto a 2 m -dimensional 
stabilizer-subspace. All stabilizer subspaces are local-Clifford 
equivalent to a linear graph-code |25|, such that 

n = c loc n G c£ c . (Di) 

Our proof utilizes this local equivalence, so first we give a 
brief account of graph codes and their relevant features. A 
graph code is defined by a graph G and a m-dimensional lin- 
ear code C over Z2. We use \G) denote the graph state corre- 
sponding to graph G, which has stabilizer generators 

kj = X s (g) Z k , (D2) 

keN(f) 

where N(j) denotes the set of vertices in the graph connected 
to vertex j. The subspace for the graph-code is spanned by 
orthogonal graph states 

\G C )=Z C \G), (D3) 

where c are binary vectors in the code C, and Z c = Zj 3 . 
The projector onto the graph code subspace is then 

U G = J2\G C )(G C \. (D4) 

cec 

For our purposes we need to express this projector in terms of 
the graph code stabilizer S 

n G = — — V s. (D5) 

seS 

The stabilizer of the graph code is 

s = { Sy = n*f ly eC± }> ( D6 > 

j 

where C 1 - is the dual of C. Allowing for local Clifford unitary 
operators, the stabilizer of the graph code is S' = C} oc SC} , 
and so the actual projector is 

II = — *— Y s'. (D7) 

s'eS' 
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In our proof we use the following fact: The stabilizers of 
S' have the same weights as the those in the locally equiv- 
alent graph code stabilizer, S. That is, if w(s) is weight of 
s (the number of nonidentity Pauli operators) then w(s') — 
w(Cio C sC' 1 o C ) = w(s), which holds because local Cliffords 
conjugate Pauli operators without changing their weight. We 
use this fact in combination with other features of graph codes. 

As for the relevant quantum state, this also has a Pauli de- 
composition. Using that a single T-magic state is 



1 / X + Y + Z 
ro = ^ ( 1 



V3 



it follows that n copies may be represented as 



5n = - y 9 (4 



geg 



V3 



«{g) 



(D8) 



(D9) 



where Q is the set of Pauli operators with positive phase. 
Hence, the projection probability is 



tr(nr »") 



= tr 



E *'-9 
E tr ( s '-3) 



1 

71 
i 

71 



(D10) 



s'GS',g&g 

The trace vanishes except when g.s' = ±1, and so 



tr(nr ® 



= E si S n ( s 



s'eS' 



1 

71 



(Dll) 



where sign(s') is ±1, matching the phase of s'. Clearly, an up- 
perbound is established when all signs are positive, and hence 



tr(nrn 



< 



E 

s'eS' 



1 

71 



(D12) 



Having arrived at an inequality purely dependent on the 
weights of S', we can use w(s') = w(s), to switch to the 
locally equivalent graph code. To determine the graph code 
weights w(s) we use the decomposition in terms of canonical 



generators expressed in Eq. ( D6 1, where every s y is identified 
with a binary vector y € C^~, 



w{s y ) 



w 



(D13) 



When multiplying generators together the Xj contributions 
can change into ±Yj, but never reduce in weight. Hence, 

w(s y )>w(y), (D14) 

where the R.H.S. is the weight, number of 1 entries, in the bit 



string y. Combining this result with Eq. ( |D12| >, we have 

tr(Hr 5 



< 



yec 1 - 



1 

71 



«(y) 



(D15) 



This inequality now depends solely on the classical linear 
code C . All such linear codes can, up-to relabeling of bits, be 
diagonalized such that the generator matrix, M, has an iden- 
tity over the first n — m elements, such that M = [l n _ m , M']. 
Dividing the bit strings into two halves y = (y',y") = 

(j/i>")J/n-m>S/i>->S/m)' then me elements of y' are fixed 
by diagonlization of the generator matrix. Furthermore, since 
w(y) = w(y') + w(y"), we can conclude 



tr(nr ® n ) < 2 r ' 



' E 

ye{o,i}"- 



1 

71 



w(y')+w(y") 



(916) 



The weights of w(y') are fixed by diagonalization, but w(y") 
depend on features of the code. We are interested in an upper- 
bound, maximized over all possible projectors n, which can 
be achieved when w(y") = 0, and so 



max [tr(nr ® ")] < 2 r ' 



E 



ye{o,i}" 



1 

71 



»(y') 
0917) 



1 



n—m 
st 



This gives an upperbound, but it is easy to verify that it is 
saturated by measuring m qubits in the computational basis 
and postselecting on "+1" outcomes. 



Appendix E 

In our magic activation protocol we saw that a(q) states 
(defined in Theorem [3]) can be a powerful resource for activa- 
tion of bound families of states. This suggests that they may 
be a powerful resource in their own right, and that the ability 
to prepare many copies of them may enable universal quantum 
computation. This is the problem we address here, although 
for a more general class of states. We consider states of the 
form 



a(q,r) = qr lfi + (1 - q 
+ r (t , + n i) , 



2r)r ,i 



(El) 



Before continuing let us reflect on some properties of these 
states. Firstly, states in this class are always separable. Sec- 
ondly, they are correlated states except for specific values 
r = ^(1 — ^fq) which give a product state. Although this 
class of states is not completely general, any state can, by a 
suitable twirling procedure (See Appendix |F]i be brought into 
this form. 

We now outline a protocol for exploiting many copies of 
these correlated states, where the maximum achievable fi- 
delity approaches 1 as r — > 0. Since the protocol consists 
in a chain of projections, we refer to the protocol as the daisy 
chain protocol 

1. Prepare n copies of the state a(q, r), with the first pair 
as qubits A & B, pair 2 as qubits C & D, and so on; 
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FIG. 5: The "phase" diagram for correlated noise resource, a(q, r), 
as described by Eq. <[ETJ. The diagram shows the region of stabilizer 
states, and the region of resources that are known to be universal 
for quantum computing. Universality is possible if the daisy chain 
protocol achieves a fidelity, see Eq. \E4\ , that exceeds the threshold 
above which the 5-qubit code can be utilized as in Ref. 1 12 1. 



2. Measure qubit A in the computational basis, and posts- 
elect on a "+1" outcome; 

3. For qubits B & C, measure XbXc and ZbZq, and 
postselect on "-1" outcomes for both; 

4. Perform the preceding step for qubits D & E, and all 
subsequent pairs; 

5. Leave the final qubit unmeasured and discard all mea- 
sured qubits. 

After step 2, qubit B is left in a product state r(/o) with fi- 
delity 



/o 



g(i-At)+r/ rt 



(E2) 



• + 9(l-At) + (l-«-2r)/ rt ' 

which is a generalization of Eq. \\5\ . 

Step 3 results in the familiar singlet projection on the sec- 
ond and third qubits (B and C). After this projection, qubit 
D is in the state t(/i), where f\ is determined by the matrix 
equation 



h 
1-/1 



q r 

r (1 - 2r - q) 



fo 
l-/o 



,(E3) 



where we use a proportionally sign as the left most vector 
must be renormalized to obtain the fidelity. This is then re- 
peated between every (2j + l) th and (2j + 2) th qubit. After 



all n — 1 singlet projections, we find that the last qubit is left 
in the state r(/ n _i), where 

n-l 

/ a r > 

oc 



fn-l 
1 - fn-l 



(l-2r-«) 



fo 
l-/o 



The limiting behavior, for large n, of this matrix equation is 
determined by the matrix eigenvalues, /xi and and eigen- 
vectors. Whenever eigenvalues have different magnitudes, the 
matrix (as n — > 00) projects onto the eigenvector E9ll with 
the largest eigenvalue. As one would expect, when a(q, r) is 
a product state the eigenvalues are identical, but in all other 
cases there is one dominant eigenvalue which determines a 
limiting fidelity 



lim /„ = 



1 1 + tan 


"1 / 


- arctan 




.2 V 



2r 



2{r + q) - 1 



(E4) 

Convergence to this fidelity is exponentially fast in the num- 
ber, n, of copies of a(q,r). Specifically, for large but finite 
n, deviations from this fidelity vanish as (/^/Vi)™, where \i<i 
is the smaller eigenvalue. In Fig. [5] we chart out various pa- 
rameter regimes indicating when the resource is a stabilizer 
state and when it provides a resource for universal quantum 
computing. Universality may be achieved by a combination 
of the daisy chain protocol followed by the standard 5-qubit 
distillation procedure lfT2l . 

Notice that product states are not the only nonstabilizer 
states where we observe a regime where no known methods 
enable universal quantum computing. There is a temptation to 
conjecture that some notion of boundness applies to any fam- 
ily of states that transverses the gap anywhere except via the 
origin. However, the daisy chain protocol shows that all non- 
product states can be purified towards some state, even if that 
state is not above the threshold for the 5-qubit code. As such 
our current definition for boundness would not extent to these 
families. However, this seems like a failing of our definition 
more than anything else, as the extent of possible purification 
appears to be limited by how correlated the raw resource is. 

Finally, note that a very similar protocol and analysis can 
be performed for 2-qubit correlated states in any basis, not 
just the T-basis. 



Appendix F 

Here we outline twirling protocols for bringing an arbitrary 
state into the form a(q, r), as defined in Eq. (El I. We perform 
the following: 



1. randomly choose a unitary operator from the set 
{1, T, T 2 } and apply to qubit A; 

2. randomly choose a unitary operator from the set 
{1, T, T 2 } and apply to qubit B; 

3. randomly choose a unitary operator from the set 
{1, Y A H A SWAP AtB H A Y A } and apply. 

The first two steps diagonalize the state in the \Tij) basis, and 
the third step mixes the symmetric terms. 
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